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ABSTRACT 


| UPDATED POSITION AND ICE VELOCITY 
FOR THE AIDJEX MANNED CAMPS 


A.S. Thorndike and T.O. Manley 


\ 


During 1975 and 1976, measurements of sea ice motion were made as 


part of the Arctic Ice Dynamics Joint Experiment (AIDJEX). In an earlier 


report (Thorndike and Cheung, 1977), energy at the inertial period was at- 
tenuated by approximately 50% due to the choice of a low covariance factor, 
Q = 100, during Kalman filtering. 


For oceanographic investigations, the data set proved to be unac~ 


ceptable particularly in the summer months when there is a significant amount 


of energy at the inertial period. 


This report is the result of the reprocessing of the original position 
data taken from the four manned camps (Caribou, Blue Fox, Snowbird and 


Big Bear) with a covariance factor Q = 100000. 


The higher Q allows data of the inertial period and longer to pass 
through the filter with less than 1 percent attenuation. 


Because of the large amount of information contained in the data set, 


it was divided into two parts. This report (Volume 2) provides position and 


ice velocity from 18 October 1975 to 4 May 1976. Volume 1,provides position 


and ice velocity from 11 April 1975 to 17 October 1975 (Thorndike and Manley, 
1980). 
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Introduction 


The objective of the Arctic Ice Dynamics Joint Experiment, as set 
forth in its most general statement, was to reach an understanding of the 
dynamic and thermodynamic interaction between the ice cover and its en- 
vironment (Maykut et al., 1972). One component of the experiment was an 
array of drifting stations at which frequent measurements of the position 
were made (Figure 1). For the present report those measurements have 


been edited using Kalman smoothing techniques to give the best estimates of 


position and velocity. (The raw measurements themselves are not reported 


here, but are available through the Polar Science Center.) The data set 
described here is restricted to the main AIDJEX experiment, which began in 
March 1975 and ended in May 1976. 

This data report represents the most complete documentation of ice 
motion available, but it is limited to one year and to the Beaufort Sea. Par- 
ties interested in earlier data, more widely dispersed in space and time, 
should contact the Polar Science Center, University of Washington, Seattle, 
Washington, which has compiled most of the ice drift data beginning with the 
drift of the FRAM in 1893 to the present. 
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Data Acquisition 


A full description of the measurement systems used in this program 
is beyond the scope of this report. The list of references includes sources 
of more detailed information (Wasilew and Vivian, 1976; Martin, P., 1978). 
All of the hardware we used was supplied by the Satellite Positioning Cor- 
poration, Houston, Texas. 

Four manned camps -- identified by radio call names Big Bear, Cari- 
bou, Blue Fox and Snowbird -- were tracked using NavSat (Navy Navigation 
Satellite System). | The geographical position of each of the four camps was 
measured at unevenly spaced times using Doppler satellite navigation tech- 
niques. The raw measurement is the frequency of a signal transmitted from 
the satellite to the ice station. The measured frequency is affected by a 
Doppler shift related to the rate of change in distance between the satellite 
and the ice station, which is itself related to the unknown ice coordinates 
and the known satellite orbit. During each satellite orbit several frequency 
measurements were made and from them an overdetermined solution was found 
for the unknown coordinates. Data recording in the field was on magnetic 

tape with hard copy printout for backup. 

A constellation of six polar orbiting satellites was maintained for this 
system. The orbital period is about 107 minutes. An observer at high lati- 
tudes can track every orbit. Of the 81 potential passes each day, many are 
unusable (because of radio interference between overlapping passes, for in- 
stance). We averaged about 25-30 usable passes per day at each station. 


Single channel (400 MHz) receivers were used to count the Doppler cycles 


for roughly 30-second intervals. The duration of each Doppler counting 


interval was measured independently. The frequency standard used in the 
receivers was stable during a 20-minute pass to about 1 part in 1012, 
Fix solutions were done in real time at each manned camp. The 


solution used only the Doppler data from the 5 minutes on either side of the 


time of closest approach of the satellite. No additional editing was performed. 


A correction to each Doppler count was based on the duration of the counting 


interval to account for the random variations in the counting interval. The 


correction also included the fixed timing delay in each receiver. The fix 


solution ignored the effect of the velocity of the ice camps. 


In fact, two receivers with antennas separated by about 100 meters 


were operated at each manned camp. In addition to serving as a backup, 


the measurements at the second antenna were used to define the azimuth of 


the line joining the two antennas. This was possible when both receivers 


produced acceptable fixes from the same satellite pass. Estimates of the 


azimuths obtained at each of the manned camps were reported by Thorndike 


and Cheung, 1977. 
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Preprocessing 


The original data tapes from the manned stations measurement sys- 
tems typically contained several "glitches" per tape because of power failures 
at the camp, or operator error, or some other reason. These ‘apes were 
transcribed discarding those records which had either a parity error or an 
incorrect number of words. The tapes which resulted from this process 
were read without error on the CDC 6400 at the University of Washington, 
where the bulk of the processing was done. 

Three problems had to be dealt with to prepare the sequence of 
measurements at each station for the actual Kalman filtering. First, data 
recording problems had created occasional gaps in the manned camp data. 
At times the backup teleprinter records from the field could be used to fill 
these gaps. This was done by keypunching the fixes directly from the 
printed records. Second, the data sequence had occasional points which 
appeared to be out of sequence in time. This was sometimes correctable, 
as when the operator had keyed in the wrong date, but usually the situation 
was ambiguous and points out of chronological order were dropped. The 
third and most serious problem involved bad fixes. A glance at a plot of 
the raw time series (Figure 2) illustrates the problem. Because outlying 
fixes have a deleterious effect on the Kalman filtering results, an attempt 
was made to remove them during preprocessing. The algorithm employed 
compares each fix in latitude with the median of the latitude of the 10 fixes 


preceding and the 10 fixes following it. When the differences exceeded a 


preset tolerance of 10 kilometers, the fix was eliminated. The algorithm 
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was applied twice, with a smaller tolerance of 1 kilometer the second time. 


The same procedure was then applied for longitude. 


Visual checks were an important part of the procedure. Plots of 


the data points were produced before and after preprocessing. At times it 


seemed from the plots that we were not detecting certain bad fixes, or that 


we were throwing out good ones, and it was necessary to adjust the toler- 


ance levels. Edited raw data plots like the one shown in Figure 3 provided 


assurance that the data were ready for filtering. 


For the main experiment, AIDJEX adopted the convention of number- 


ing days consecutively, beginning with day 1 = 1 January 1975 and ending 


with day 500 = 14 May 1976. Throughout this data report, this numbering 


system is used quite extensively. Appendix 3 provides a conversion table 


from AIDJEX days to calendar days. 
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Kalman Smoothing 


In the smoothing used here an assumption is made regarding the 
motion of the ice. If the state of the ice at time t, is represented by a 
vector Xn , containing, say, the position, velocity, and acceleration of a 
point, then the state at some future time is assumed to be partially deter- 


mined by Xn : 


Xnet = OX, +TWoay (1) 


where ¢ and [I are matricies known from the physics of the situation, and 
Wh +1 represents an unknown random perturbation. In Kalman's formulation, 


which we follow, the random effects are assumed to be Gaussian and White: 


Q ifn=m, and 
cov (W_W_) = (2) 
me 0Oifn#m 


2, 76 


We have chosen Q = 100000 m“hr . Qualitatively choosing a small value of 
Q is equivalent to assuming that the ice moves in a completely predetermined 
way. (See from equations 1 and 2 that Q = 0 would imply that xy Po gs oX . 
a large value of Q corresponds to assuming that the ice experiences large 
sudden changes in its acceleration.) Thorndike and Cheung (1977) used 

Q = 100 m?hr®, The lower value of Q assumed smooth behavior, but was 
unnecessarily conservative. The present value of Q was determined by trial 
and error. A sequence of raw fixes was used as input to the Kalman filter. 


Several different runs were done using different values of Q. The present 


value of Q appeared to give the best agreement between the output from the 


filter and the raw fixes. 


The original Q (100 mhr~*) used by Thorndike and Cheung removed 
approximately 50% of the inertial period amplitude as can be seen in Figure 4. | 
A constant amplitude sinusoidal function whose period varied linearly with 


time was chosen as the input to the Kalman filter with Q = 100 m?hr®, In- 


put periods ranged from 0.5 hours (at day 10) to 60.4 hours (at day 15) with 


the inertial period being located at day 10.9. An approximate 50% attenuation 


of the input function at the inertial period can be roughly estimated from the 
graph. 


Figure 5 shows the same input function and the resulting Kalman 


output with Q chosen to be 100000 m hr-®, The inertial period passed through 


the filter with very little attenuation. Another test was then run, this time a 
constant amplitude, constant period (12 hours) sine wave was used to simulate 


the inertial period of the Arctic latitudes. The results indicated that there 


was only a 0.03% attenuation from the original input function. 


For the Kalman reduction, we defined A = th ae? t 


n? and the 
matrices 
x = position @ = #1 =A a2) 2 
velocity 0 1 A (3) 
acceleration 0 0 0 
r= |a%e6 
7/2 (4) 
A 
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The dimensions of the matrices are X : 3x1, Q:3x3,T:3x1, 6: 3x 3, 
W:1x1, andQ:1xl. 

An identical equation can be written for the y-coordinates. The 
two coordinates are treated independently. 


The measurement Z = is related to X by 


Zn = (1 0 Ox, +V, 


Z 3 : 2 
where vA is a measurement error with assumed covariance R = (60 meters) . 


The response of the measurement and processing system was investi- 
gated as follows. Four hundred sample points were drawn independently 
from a Gaussian distribution with zero mean and variance (50 meters)”, A 
sampling interval of one hour was assumed. This noise simulated the Nav- 
Sat measurements. The sequence was used as input to the Kalman filter 
procedure. The velocity output from the filter shows two effects. First, 
since velocity is the time derivative of position, we expect to find the spec- 
trum S,(%) behaving as ws ( w), where S,(w) is the spectrum of position. 
S$, (%) was assumed to be white; 8, () = o” , which implies Ss, () = we a 5 
The second effect is the tendency of the filter to smooth out high frequen- 
cies. Qualitatively then, we have something like the diagram in Figure 6. 

The actual spectrum of ice velocity is shown in Figure 7. The dia- 
gram illustrates how the measurement errors and finite sampling rate limit 
the accuracy of the velocity estimates. Most of the variance in the velocity 
errors occurs for the frequencies greater than 3 cycles per day. For lower 


frequencies we can estimate the contributions to the variance as 0 - 1.5 cycles 
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per day: o2 = 3x 108m? sec 2, 1.5- 3.0 cycles per day: ot 2x 107° 


m” sec™” total variance 0 - 12 cycles per day: o” = 104m? sec 2. 

The general smoothing problem is to estimate X at time th using the 
entire set of observations (Z;) i=1....m, and to give the error variance 
of that estimate. The solution equations, due to Kalman, are reproduced 
in Appendix 1. All positions were converted to an X-Y coordinate system 
prior to actual filtering. This coordinate system along with its governing 


equations are shown in Appendix 2. Typical estimation error variances for 


position and «--xcity are 40 meters and 1.0 cm/sec respectively. 
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Results 


The final data listing consists of five major columns. The first 
column is the time in decimal AIDJEX days corresponding to the observation 
of the position and ice velocity at a particular camp. The remaining columns 
indicate the camps for which data are applicable. Subdivisions under each 
major camp heading are LAT (latitude in decimal degrees) , LONG (longitude 
in decimal degrees), and the north and east components of ice velocity in 
Cm/sec, NIV and EIV respectively. 

Due to the large amount of information contained in the data set, it 
was divided into two parts. This report (Volume 2) covers AIDJEX days 
291 (18 October 1975) to day 490 (4 May 1976). Volume 1 covers AIDJEX 
days 101 (11 April 1975) to 290 (17 October 1975) (Thorndike and Manley, 
1980). 

Figures 8 - 11 show a more detailed drift track of each of the manned 
camps. An insert shows where the region lies with respect to the north slope 
of Alaska and Banks Island. Asterisks superimposed on the drift track indi- 
cate the positions at integral multiples of 20 days. The beginning and end- 
ing days are noted for each trajectory. A dashed line indicates a long period 
of missing data. The plotted region is 500 x 500 km and is aligned with the x, 
y coordinate system shown in Appendix 2. Figures 8 - 11 were reproduced 
from Thorndike and Cheung (1977). 

Table 1 lists the last known positions and ice velocities from camp 


Big Bear during the spring of 1976. Big Bear was the original main camp 


but was abandoned in early October of 1975 due to ice breakup. In spring 


of 1976, RAMS buoy number 604 was placed on the site of the camp which by 
this time had been spread over a few kilometers. Data cover from day 465 


to 485. Because of RAMS buoy lower sampling rate and higher noise prob- 


lem, the data in Table 1 were not reprocessed. 
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APPENDIX 1 


KALMAN SMOOTHING EQUATIONS 


state vector (n x 1) 
transition matrix (n xn) 


coefficient matrix for random perturbations tc physical 
model (n x m) 


white noise driving the physical process (m x 1) 
covariance matrix for W (mx m) 

measurement vector (2 x 1) 

matrix relating measurement to state vector (2 xn) 
measurement error vector (2 x 1) 

measurement error covariance (2x 2) 


prediction, filtering and smoothing error covariance 
matrix (n x n) 


identity matrix (n xn) 

Kalman gain matrix for filtering (n x £) 
smoothing gain matrix (n x n) 

number of measurements in the data sequence 
index for the data sequence 

time 


The superscript T denotes the transpose operation. The symbol 


J 


The notation here follows Meditch (1969). 


Prediction 


Given observations ( 2 sees 2)» predict the state at the future 


time tk +1 and find the prediction error covariance. 
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x(k+1,[k) = O(k +1, k) X(kK]k), A(1) 
P(k+1]k) = @(k+1,k)P(k[k) oO? (k + 1,k) 

+ I (k+1,k)Q(k) I? (k =1,k). A(2) 
Filtering 


Given observations (2), sere 2) and a new observation ZK+1) , 
estimate the present state at time te +1 and find the filtering error co- 
variance. 


XCk+1i[k+) = Le +1[k) + K(k + 1 (Zk + 1) | 
- H(k +1) X(k+1/k)), A(3) 
P(k+i]k+1) = (I1- K(k + 1)H(k +1) P(k+1]k), A(4) 
where 
K(k +1) = P(k +1] k) H! Ck + 1) (H(i + 1) P(e + 1]k) H (k + 1) 
+ Rk +1))2 ACZ 


Smoothing 


Given the complete data set (Z;, sees Zn)» estimate the state at 
any point te and find the smoothing error covariance. 


KCkIN) = XC [k) +A Ck) (XCk + 1[N) - X(k + 1/k)) “A(6) 
P(k|]N) = P(k|k +A(k) (P(k + 1]N) 

- P(k +1/k)) A? (k), where A(T) 

A(k) = P(k[k) ¢7(k + 1, k) Plc + 1] k). A(8) 
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Implementation of the Filter Algorithm 

Starting from an initial condition x (0| 0), P(0| 0), a prediction is 
made to time th: The prediction error covariance matrix P (1 | 0) is used to 
obtain the Kalman gain K(1). The filtered estimate X(1| 1) is then found 
by adding to the predicted x (1| 0) a quantity which is proportional to the 
difference between the observation Z(1) and H(1)X(1]0). H(1)X(1[0) is 
the value the measurement would have if the true state were actually X(1 | 0) 
and no measurement errors were made. The algorithm is repeated as each 
new data point is entered until the last filtered quantities X (N|N) and 
P(N |N) are produced. 

The smoothing algorithm, proceeding backward in time, utilizes 
results from the prediction and filtering calculations. The first step is to 
find X (N-1|N) and P(N-1|N) from equations A(6), A(7) and A(8). No 
direct use of the measurements is made in smoothing. The algorithm simply 
corrects the filtered vector x (k|k) by a quantity proportional to the differ- 
ence between the prediction vector at time tk +1 and the smoothed value at 
time tet’ 

In the calculations done for this application, the CDC 6400 computer 
time required to do the prediction, filtering, and smoothing is about 0.5 sec- 
ond per data point, nearly half of which is absorbed in the matrix inversion 
in A(8). Because of the recursive nature of the filter, demands on central 


memory storage space are very modest. 


Reproduced from Thorndike, A.S., 1973, An integrated system for meas~ 
uring sea ice motions. In Ocean 73, IEEE International Conference on En- 
gineering in the Ocean Environments, IEEE publ. no. 73 CHO 774-0 OEC, 
pp. 490-499. 
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APPENDIX 2 
COORDINATE SYSTEMS 


Position measurements were made in geographical coordinates 
(latitude north, longitude east). The smoothing operation was 
done in a Cartesian system (X,y), where 


x 
y 


110.949 (90° - latitude) cos Qongitude) (km) 
110.949 (90° - latitude) sin (longitude) (km) 
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APPENDIX 3 


CONVERSION TABLE FOR AIDJEX DAYS 
TO CALENDAR DAYS 


For the main experiment, AIDJEX adopted a convention of number- 
ing days consecutively, beginning with day 1 = 01 January, 1975 and ending 
with day 500 = 14 May, 1976. 

In the conversion table, the first column is the AIDJEX day, the 
second is the corresponding day of 1975 or 1976 and the third entry is the 
calendar date. 


Appendix 3 - Conversion Table 
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TABLE 1 


Last known positions and ice velocities for camp Big Bear. Taken 


from RAMS buoy R604 (Thorndike and Cheung, 1977). 
AIDJEX days. 


Days are in decimal 
Negative longitudes imply west longitude. NIV and EIV are 


north and east ice velocities respectively. Negative values in the NIV or 


EIV columns imply south or west velocities. 
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Day 
465.50 
465.75 
466.00 
466.25 
466.50 
466.75 


467.00 
467.25 
467.50 
467.75 
468.00 
468.25 
468.50 
468.75 


469.00 
469.25 
469.50 
469.75 
470.00 
470.25 
470.50 
470.75 


TABLE 1 


Last known positions and ice velocities for camp Big Bear. 


Taken from RAMS buoy R604. 
(Thorndike and Cheung, 1977) 


Latitude 


72.9569 
72,9563 
72.9552 
72.9496 
72.9394 
72.9276 


72.9167 
72.9089 
72.9057 
72.9068 
72.9121 
72.9214 
72.9272 
72.9299 


72.9297 
72.9286 
72.9302 
72.9327 
72.9336 
72.9326 
72.9306 
72,9278 


Longitude 
-141.1691 


-141. 1278 
~141.1166 
-141.1382 
-141.1717 
-141.1990 


~141.2236 
-141.2530 
-141.2897 
-141. 3366 
-141.3963 
-141. 4686 
-141.5281 
-141.5752 


-141.6059 
~141.6105 
~141.5871 
~141.5560 
-141.5304 
-141.5031 
-141. 4635 
-141.4068 
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NIV EI 
ee | 8.0 
1 4.6 
~1.4 -1.2 
~4.4 -5.5 
-5.9 -4.7 
-6.2 -3.6 
-5.0 -4.0 
~2.9 ~5.0 
~ 4 ~6.3 
2.4 -8.9 
5.1 -11.5 
5.3 ~11.5 
2.9 -8.8 
0.7 -6.4 
- .7 ~2.9 
- 1 1.6 
1.6 5.5 
1.0 4.3 
- 2 3.4 
- .7 5.2 
-1.0 8.3 
-1.8 9.5 
continued - 
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Day 
471.00 
471.25 
471.50 
471.75 
472.00 
472.25 
472.50 
472.75 


473.00 
473.25 
473.50 
473.75 
474.00 
474.25 
474.50 
474.75 


475.00 
475.25 
475.50 
475.75 
476 .00 
476.25 
476.50 
476.75 


Latitude 


72.9226 
72.9166 
72.9143 
72.9180 
72.9293 
72.9385 
72.9459 
72.9513 


72.9550 
72.9580 
72.9619 
72.9646 
72.9640 
72.9615 
72.9587 
72.9558 


72.9521 
72.9459 
72.9404 
72.9370 
72.9326 
72.9251 
72.9182 
72.9188 


TABLE 1 (cont.) 


Longitude 
-141.3516 
-141.2966 
~141.2500 
-141.2199 
~141.2116 
-141.1965 
-141.1772 
-141.1618 


-141.1577 
-141.1604 
~141.1618 
-141.1677 
-141.1849 
-141.2068 
~141.2205 
-141.2203 


-141.2160 
-141.2160 
-141.2185 
-141.2307 
-141. 2526 
-141.2885 
-141.3358 
~141. 3874 
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NIV EIV 
-3.2 8.8 
-2.5 8.1 
.2 6.1 
5.3 1.7 
6.9 4 
4.5 2.5 
3.1 3.1 
2.4 1.6 
1.3 - .5 
1.7 - .6 
2.3 1 
5 -1.9 
-1.0 -3.3 
-1.5 ~3.1 
-1.4 - .8 
-1.6 1.0 
-3.1 7 
-3.7 0.0 
-2.6 -1.1 
-1.2 -2.3 
-3.7 -4.3 
-6.8 -6.8 
-2.1 -7.7 
2.8 -8.2 
continued - 


J 
¢ 
TABLE 1 (cont.) 
Day Latitude Longitude NIV EIV 
477.00 72.9276 -141. 4327 7.0 -5.7 
477 25 72.9388 -141. 4558 4.4 -1.4 
477.50 72.9441 -141.4542 Ll 1.4 
477.75 72.9443 ~141.4557 -1.3 ~1.5 
478.00 72.9406 -141.4666 -2.5 -1.8 
478.25 72.9345 -141.4747 -3.9 - .5 
478.50 72.9255 -141.4718 -5.4 1.4 
478.75 72.9170 -141.4579 -3.6 3.0 
479.00 72.9123 -141.4397 -1.1 2.7 
479.25 72.9125 -141.4319 2.4 -1.3 
479.50 72.9174 -141.4436 2.7 -2.3 
479.75 72.9195 -141.4505 .4 - 7 
480.00 72.9192 -141.4565 ~ .8 -1.0 
480. 25 72.9183 ~141.4761 ~ 1 -5.0 
; 480.50 72.9191 -141.5241 9 -9.5 
z 480.75 72.9218 -141.5924 1.9 -11.4 
481.00 72. 9266 ~141.6676 3.1 -11.1 
q 481.25 72.9324 -141.7305 3.1 ~4.9 
4 481.50 72.9393 ~141.7560 4.1 -3.3 
5 481.75 72.9489 -141.7819 5.8 -4.5 
4 482.00 72.9674 -141.8250 8.0 -8.6 
F 482.25 72.9790 -141.8889 8.9 -10.6 
4 482.50 72.9965 -141. 9578 8.9 -10.5 
. 482.75 73.0138 -142.0193 8.9 -8.2 
i ' continued - 
| 
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Day 
483.00 
483.25 
483.50 
483.75 
484.00 
484.25 
484.50 
484.75 


485.00 
485.25 


Latitude 
73. 
73. 
0496 
73. 
73. 
73. 
73. 
0542 


73 


73 


73. 
.0373 


73 


0310 
0449 


0499 
0506 
0543 
0561 


0487 


TABLE 1 (cont.) 


Longitude 
-142.0579 
-142.0749 
-142.0779 
-142.0832 
~142.0979 
-142.1235 
-142.1321 
-142.1247 


~142.0968 
-142.0555 
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Position and Ice Velocity Data 


_ 
. The following section comprises the updated positions and ice veloci- 
; ties for the four manned camps Caribou, Blue Fox, Snowbird and Big Bear. 
4 Data span from 18 October 1975 to 4 May 1976. 
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